Dynamics of a classical gas including dissipative and mean field effects 
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By means of a scaling ansatz, we investigate an approxi- 
mated solution of the Boltzmann-Vlasov equation for a clas- 
sical gas. Within this framework, we derive the frequencies 
and the damping of the collective oscillations of a harmoni- 
cally trapped gas and we investigate its expansion after release 
of the trap. The method is well suited to studying the colli- 
sional effects taking place in the system and in particular to 
discussing the crossover between the hydrodynamic and the 
collisionless regimes. An explicit link between the relaxation 
times relevant for the damping of the collective oscillations 
and for the expansion is established. 



The favored signature of Bose-Einstein condensation 
in weakly interacting gases is the time-of-flight expansion 
[1]. In this technique, the asymmetric trapping potential 
is switched off and the evolution of the spatial density 
is monitored. After long time expansion, the observed 
inversion of the aspect ratio reflects the anisotropy of the 
initial confinement. In an ideal Bose-Einstein condensate 
(BEC), this effect is a direct consequence of the Heisen- 
berg uncertainty constraint on the condensate wave func- 
tion. For an interacting BEC, the inversion is also pro- 
duced by the anisotropy of the pressure gradients caused 
by the hydrodynamic forces. The changes in the shape of 
the expanding gas can be characterized by scaling factors, 
which provide an easy quantitative tool for the analysis 
of on-going BEC experiments. The set of equations for 
those factors have been derived in many papers [2,3]. A 
similar effect has also been predicted for a Fermi gas in 
its superfluid phase [4]. Strong anisotropy has been re- 
cently measured in the expansion of a highly degenerate 
Fermi gas [5] close to a Feshbach resonance. Resonance 
scattering can also give rise to anisotropic expansion in 
the normal phase as proven in the experiments of [6,7] 
carried out in a less degenerate regime. Some recent ex- 
periments on bosonic atoms above the critical tempera- 
ture have also reached the collisional regime investigating 
both the oscillations of the low lying quadrupole mode 
and the expansion in asymmetric traps [8,9]. 

So far analytic calculations for the expansion of a clas- 
sical gas have been limited to either the ballistic or to the 
hydrodynamic regime [3] . As a consequence, it is impor- 
tant to generalize such calculations to all intermediate 
collisional regimes. This is precisely the main purpose 
of this paper. We begin by an outline of the theoretical 
description of the thermal gas based on the Boltzmann- 
Vlasov equation. Our approach relies on an approxi- 



mated solution to this equation by means of a scaling 
ansatz. This solution is used throughout the paper to 
investigate two kinds of related problems: the lowest col- 
lective oscillation modes and the time-of-flight expansion 
when the confinement is released. 

The Boltzmann-Vlasov (BV) kinetic equation for the 
phase space distribution /(t, r, v) takes the form: 
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where Uh is the trapping potential chosen of harmonic 
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Interparticle interactions en- 



ter Eq. (1) in two different ways [10]. On the one hand, 
they modify the effective potential through the mean field 
term U m [ which affects the streaming part of the Boltz- 
mann kinetic equation. The mean field potential U m { is 
equal to 2gn for bosons and gn for two fermion species 
[11], where the coupling constant g — iirh 2 a/m is fixed 
by the s-wave scattering length a. The mean field term is 
linear in a and is non dissipative. On the other hand, two 
body interaction determines the collision integral / co n [/] 
which is quadratic in the scattering length, and describes 
dissipative processes. Eq. (1) is valid in the semi-classical 
limit, namely when the thermal energy is large compared 
with the separation between the energy eigenvalues of the 
potential [12,13]. 

In this paper, we will treat the collision integral within 
the relaxation time approximation. This model should 
suffice to capture the essential physics of the problem. 
We consequently write: 
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where r is the relaxation time related to the average time 
between collisions, and f\ c is the local equilibrium den- 
sity in phase space. As a consequence, f\ e has a spherical 
symmetry in velocity space, i.e. it depends on the veloc- 
ity through [v — u(r)] 2 where u(r) is the local velocity 
field. 

The dynamics of the gas will be described by the fol- 
lowing scaling ansatz for the non equilibrium distribution 
function: 
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where /o is the equilibrium distribution function which 
satisfies the equation (7 co ii[/o] = 0): 
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The dependence on time of / is contained in the dimen- 
sionless scaling parameters bi and 9i. The parameter 6, 
gives the dilatation along the i t;i direction, while #i gives 
the effective temperature in the same direction. Such an 
ansatz generalizes the one used in [14]. We recall that 
in this method the shape of the cloud does not explicitly 
enter the equations. This is the reason why the solutions 
are equally valid for a dilute Bose gas above the critical 
temperature, a dilute Fermi gas in its normal phase, and 
a classical gas. 

Following Ref. [14] , one can derive the set of equations 
for the scaling parameters bi and 9% (see Appendix A): 
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where the dimensionless parameter £ = (J7 m f)o/((^mf)o + 
2m{v 2 )o/3) accounts for the mean field interaction [15] 
and 9 = ^ . 0j/3 is the average temperature, for a clas- 
sical gas (v 2 )o = SksT/m. The parameter £ is expected 
to be small for dilute gases (na 3 c 1) since the ratio 
U m t/k B T scales as (na 3 ) 1 / 3 (nA^ b ) 2 / 3 , where \ db is the 
de Broglie wavelength and n the mean density [14]. Eq. 
(6) shows that the dissipation occurs when the tempera- 
ture is not isotropic and the relaxation time r has a finite 
value. 

Eqs (5) and (6) are the main result of this paper. The 
collisionless regime is obtained by taking To = oo. In this 
limit, we have the simple relation 9- L = b~ 2 between the 
scaling parameters, and we recover the equations derived 
in [14]. In the opposite limit (hydrodynamic regime), lo- 
cal equilibrium is always ensured because of the high col- 
lision rate. As a consequence, we have 9i = 9 = . b^ 
and the Eqs (5) and (6) can be recast in the form: 
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For £ = (no mean field), we recover the equations first 
derived in [3]. Note that in both the collisionless and 
the hydrodynamic regimes, the collisional term does not 
contribute since there is no dissipation in these limits. 
We next focus our attention on the intermediate regimes 
where the collision term enters explicitly the equations 
of motion. 

Let us first study the breathing mode in the case of 
spherical harmonic trapping with angular frequency u>o- 
In this case, we find a solution with bi = b and 9i = b~ 2 . 
For such a solution the collision term identically vanishes 



in all intermediate collisional regimes. Our approach can 
be readily generalized to lower dimensions leading to the 
frequency ojo(4 + £(d — 2)) 1 / 2 for the monopole mode 
[14], where d is the dimension of space. In two dimen- 
sions the mean field does not affect the frequency of the 
monopole. This comes out from the fact that in this case 
the ansatz is an exact solution of the BV equations, as 
already stressed in Ref. [16]. 

We now consider a sample of atoms confined in a three- 
dimensional cylindrically symmetric harmonic potential. 
We denote by A = uo z /uu^_ the ratio between the axial 
and radial angular frequencies. Expanding Eqs. (5) and 
(6) around equilibrium (bi = 9i = 1) we get a linear and 
closed set of equations which can be solved by looking for 
solutions of the type e %u>t . The associated determinant 
yields the dispersion law: 
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where A[lu] = lu 2 (lu 2 — lu 2 1+ )(uj 2 - w 2 j_), B[lu] = uj(lu 2 — 
^hd+X^ 2 - w hdJ, CM = u{uj 2 - uj 2 c1 ) and D[lu] = 
(u 2 — w^ d ) and To is the value of the relaxation time 
r calculated at equilibrium and 



c4 ± = ^4(l + A 2 )-A 2 e 

±v/16 + A 4 (4 - i) 2 + 8A 2 (£ 2 - 4 + o) 
a; 2 =a,i(4-20 
<i ± = ^ (5 + 4A 2 + £(1 + A 2 /2) 

±^(10 + 8A 2 + 2£ + A 2 £) 2 - 72A 2 (4 + f ) 

Here £ is the parameter accounting for the mean field 
effects and (cl) and (hd) refer to the collisionless and 
hydrodynamics regimes respectively. The solution of Eq. 
(8) interpolates the frequencies of the low lying modes 
for all collisional regimes ranging from the collisionless 
to the hydrodynamic. As the confinement is cylindrically 
symmetric around the z axis, we can label the modes by 
their angular azimuthal number M. The first factor of 
the l.h.s. of Eq. (8) gives the frequencies of the two 
M = modes, while the second factor gives that of the 
quadrupole (M = ±2) mode. The roots of A and C 
have already been obtained in [14], and correspond to 
the frequencies of the low lying modes of a collisionless 
gas in presence of mean field. Eq. (8) for £ = has been 
derived in Ref. [17] and the corresponding frequencies 
have been investigated experimentally [8]. For £ = 1, 
corresponding to (gn )o > (v 2 )o, we find lu 2 1± = w 2 d± , 
cj 2 i = w^ d , and the frequencies coincide with the ones 
predicted for a Bose-Einstcin condensate in the Thomas- 
Fermi regime [18]. 
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So far, we have not given the explicit link between the 
relaxation time entering Eq. (8) and the collision rate. 
Following Ref. [17], we can establish this link for a classi- 
cal gas by means of a Gaussian ansatz for the equilibrium 
distribution function /(r,v, t). One obtains tq = 5/(47) 
where 7 = 2(27r)~ 1 / 2 n max ert;th is the classical collision 
rate where vth = (ksT /m) 1 / 2 is the thermal velocity, 
"max is the peak density and a is the cross section which 
is assumed to be velocity independent. For bosons the 
link between the scattering length and the cross section 
is a = 8ira 2 whereas for two fermions species one has 
a = Ana 2 . 

We now establish the set of equations that describe the 
time-of-flight expansion. In the collisionless regime where 
the mean free path is very large with respect to the size of 
the trapped cloud and in the absence of mean field contri- 
bution, we readily obtain the exact equations hi = ujf /bf 
which admit the solutions 6j(t) = (1 + ujft 2 ) 1 / 2 , leading 
to isotropic density and velocity distributions after long 
time expansion. 

When the effect of collisions is important the physics 
of the expansion changes dramatically. As an exam- 
ple, the radial directions of a cigar-shaped cloud expand 
faster than the longitudinal one finally resulting in an 
anisotropic velocity distribution. So far, an analytic ap- 
proach has been proposed only in the full hydrodynamic 
regime [3]. However, this approach assumes that the hy- 
drodynamic equations are always valid during the expan- 
sion. In general, this cannot be the case since the den- 
sity decreases during the expansion reducing the effect of 
collisions. Alternatively, the expansion of an interacting 
Bose above T c gas has been investigated by means of a 
Monte Carlo simulation [19]. 

In our approach, we provide an interpolation between 
the two opposite collisionless and hydrodynamic regimes 
using the scaling formalism. The decrease of the collision 
rate during the expansion yields a non constant relax- 
ation time r(6j, 6{) that depends explicitly on the scaling 
parameters and reflects the changes of the density and 
the temperature during the expansion. As a result, the 
expansion is described by the following set of 6 non- linear 
equations 
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The dependence of the relaxation time r on the scaling 
parameters is obtained by noticing that the collision rate 
7 scales as nT 1 ! 2 . Using the scaling transformation n — > 
n o(nj bj) and T — > Tq9, where no and To are the initial 
density and temperature respectively, we deduce 
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where tq is the average time of collision at equilibrium 
[20]. Since both the results (8) for the dispersion of the 
linear oscillations and Eqs. (9) and (10) for the expansion 
have been derived starting from the same scaling equa- 
tions (5) and (6) the relaxation time To entering the two 
processes is the same. As a consequence, the combined 
investigation of the expansion and of the quadrupole os- 
cillations can provide a useful check of the consistency 
of the approach and, possibly, useful constraints on the 
value of the cross section. 

The time evolution of the aspect ratio R± (t) / R z (t) = 
^b±(t)/b z (t) in the absence of mean field is depicted on 
Fig. 1 for different values of the product lu±tq and for 
an initial aspect ratio A = 0.1. In the collisionless regime 
(to = 00), the aspect ratio tends asymptotically to unity 
reflecting the isotropy of the initial velocity distribution. 
For other collisional regimes, the asymptotic aspect ratio 
is larger that one as a consequence of collisions during the 
expansion. We find a continuous transition from the col- 
lisionless to the hydrodynamic prediction as we decrease 
To from infinity to zero. We then conclude that in general 
the expansion cannot be described with either the hydro- 
dynamic or the collisionless prediction [21], but requires 
a full solution of our equations (9). Similar conditions 
have been already encountered experimentally for classi- 
cal or almost classical gases [9] . We also notice that it is 
very important to take into account the time dependence 
of the relaxation time, accounted for by the scaling law 
(10). For example by simply using t = tq during the 
whole expansion, the curve lu±t = 0.1 of Fig. 1 (solid 
line) would be shifted upward and the resulting predic- 
tion would result much closer to the hydrodynamic curve 
(dotted line). 
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1. Aspect ratio as a function of the normalized 
time u±t for different collisional regimes (initial aspect ratio 
A = 0.1): collisionless (dashed line, to — > 00), intermedi- 
ate collisional regime (solid line, u±tq — 0.1 J, hydrodynamic 
regime (dotted line, tq — > 0). 



Let us finally comment on the effect of quantum statis- 
tics on the calculation of the relaxation time. For a Bose 
gas at temperature above T c the problem has been inves- 
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tigated in [22] where it has been shown that statistical 
effects do not play a significant role. In contrast, the 
relaxation time in a harmonically trapped dilute Fermi 
gases has been shown to be strongly affected by Pauli 
blocking at low temperature [23]. The effects of colli- 
sions in a strongly interacting Fermi gas, including the 
unitarity limit, have been recently addressed in [24]. 



obtain a relation among the 6 scaling parameters one 
uses the identity (v 2 ) — (v 2 )\ c , from which we deduce 
6 = The contribution of the local equilib- 

rium term to the second term is obtained by noticing 
that, at lo_cal_ equilibrium, 9f _= 6: J V?f ie d 3 Rd 3 V = 
T / V 2 f (R,V)d 3 Rd 3 V = NT9k B T /{m9i). Hence Eq. 
(A2) can be recast in the form (6). 
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APPENDIX A: EQUATIONS 

We make the following ansatz for the non equilibrium 
distribution function: /(r,v,i) = r/o(R(t), V(t)) with 

Ri = n/b t , v t = ( Vi - WW 1/2 and r = U 3 b J% 1/2 - 

The dependence on time is contained through the param- 
eters bi and 0j. Following [14], we substitute this ansatz 
into Eq. (4) and use the equation for the equilibrium 
distribution /q. Wc find 



1 



1 



b,e 



V2 IT A 



Ri 

^72 



bi + w, bi - 



bi Ujbj 



'coll 



(Al) 



Performing integra- 
tion in phase space, we calculate the average moment 
of RiV z , namely / RM[Eq.(Al)]<PRtjPV/N. This leads 
to Eq. (5). Note that this equation is not affected by the 
collision integral since the quantity RiVi is conserved by 
collisions. 

To derive Eq. (6), we consider instead the average 
moment of V 2 . This yields: 



NTk B T 
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(A2) 



where T is the equilibrium temperature. Differently 
from (5), Eq. (A2) depends explicitly on the collision 
integral . In order to calculate the r.h.s of Eq. (A2), we 
use the relaxation time approximation: I = — (/— /i e )/ro. 
The first term gives: / V 2 fd 3 Rd 3 V = NTk B T /m. To 
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